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1. Introduction 

Today we confront the challenge in cosmology of explaining the observation that, in the 
recent past, our Universe has experimented a phase of accelerated expansion [1]. This 
acceleration is attributed to the mysterious dark energy (for a review see pj), whose 
understanding is nowadays a major problem. Although the simplest candidate to play 
the role of dark energy is the cosmological constant A allowed by Einstein equations, its 
smallness Kobserved ~ lO^^^^Mp is very problematic as it requires an extreme fine-tuning. 
Besides, recent observations indicate that a time dependent equation of state is possible. 
In this case, the simplest approach is to consider a scalar field to play the role of dark 
energy. There are many types of such models and an extended literature, we refer only 
a few works in four of these approaches namely quintessence [3] , k-essence [1] , tachyons 
[5J and quintom [6]. 

Moreover, the SNIa data admits an equation of state Ude < — 1, which is attributed 
to the so-called phantom dark energy [7j. This has the striking feature that its density 
grows without limit with the expansion of the Universe. Usually, this behavior leads to 
a violation of the weak energy condition and then, to the so-called Big Rip in a finite 
time [8j. This is true if dark energy satisfies a conservation equation which corresponds 
to a non-interacting fluid. Now, if dark energy interacts, for example with dark matter, 
then the energy conservation equation is modified, and it is possible to circumvent the 
blowing of the dark energy density. Further, there are proposals that can encode the 
crossing ujde = — 1 (dubbed phantom divide) without violating the weak energy condition 
[9]. This behavior can be exhibited in dark energy-dark matter unified models. An 
example of such a model is the Chaplygin gas [10] and its generalization [11]. The 
Chaplygin gas has a connection with string theory and can be obtained from the light- 
cone parametrization of the Nambu-Goto action, associated to a Z)-brane [12]. In this 
model, a single self-interacting scalar field is responsible for both dark energy and dark 
matter, giving also the observed accelerated expansion. In a recent work [13], it has 
been claimed that the galaxy cluster Abell A586 exhibits evidence that dark energy and 
dark matter are interacting. The authors trace back the coupling to a departure of the 
equilibrium settled by the virial theorem and show that it could be explained by the 
generalized Chaplygin gas (GCG) with dark energy equation of state Ude = — 1- 

In the light of these last results and the considerations mentioned above, in this 
paper we will consider the generalized Chaplygin gas as a source of dark matter and 
dark energy which are interacting assuming a time dependent dark energy equation of 
state with the property to encode the phantom-like behavior. For simplicity, we will not 
consider the baryonic or other kind of matter. A particular and simple ansatz on the 
separation of the dark energy and dark matter densities is assumed. We find a phantom- 
like equation of state for dark energy that crosses the phantom divide uJde = — 1 in recent 
epoch. Its associated density smoothly decays with the expansion of the Universe (thus, 
phantom-like). It is also shown that the normalized densities fij are in the past in 
agreement with the ACDM model, but in the future they enter in a phase of mutual 
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equilibrium, where their final values are practically the same as today. This indicate 
a scaling solution, similar to that obtained in [H], where it was assumed a stationary 
dark energy which could solve the coincidence problem. Also, we find a coupling that 
it is in agreement with one would expect in relation with the coincidence problem. 

We organized the present paper as follows: In section 2, we briefly revisit the unified 
generalized Chaplygin gas. In Section 3, we construct the model by splitting the total 
density into dark energy and dark matter components, by means of our assumption. In 
section 4, we find a dark energy equation of state and we analyze its general behavior. 
We obtain also the involved densities explicitly in function of the scale factor, and show 
that the model has a scaling solution. At the end of this section, we plotted the relevant 
quantities for best fitted values of the parameters. Finally, we present the conclusions 
in section 5. 



2. The generalized Chaplygin gas 

In the GCG model the pressure p of the fiuid has the following form [TT] 

P= > (1 

where p is the total density and A > and a > are parameters. For a = 1 we recover 
the standard Chaplygin gas and a = corresponds to the ACDM model. Following the 
definition p = up, we read from ([1]) the GCG equation of state (EoS) 

u^ = -Asr-] , (2) 



p 

where Ag = A/pJ^" and po is the total density today, given by po = (3/87rG)ifo in 
terms of the Hubble constant Ho in a fiat universe. In the ranges of the values of 
As and a have been analyzed, and it is argued that the observations favor < Ag < 1 
with a > 1, although < As < 1 with < a < 1 is still possible. In this work, we will 
apply our results in both regions. The density satisfies the conservation equation 

cip + 3— (p + p) = 0. (3) 

a 

Using ([!]) in this last equation, we find the density in terms of the scale factor 



P 



1 

l + Q, 



Po, (4) 



where we normalized the scale factor as Oq = 1 today. The GCG EoS ([2]) in terms of 
the scale factor is 

= " 1 - As + A,a3(i+") • 
It is easy to see that with the restriction Ag > 0, this EoS is constrained to the interval 
—l<uj<0. Therefore, the model is in general of quintessence and excludes the 
phantom region u < —1, with a phase of acceleration for u < —1/3. Furthermore, it 
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allows us to interpret the fluid as cold matter ~ for a ^ 0, and as dark energy 
u; ~ — 1 for a ^ oo and it can be considered as an unified model of matter and dark 
energy. If we assume in the simplest case a EoS for dark energy ujde = — 1, the resulting 
model gives slight deviations of the ACDM model [11]. However, due to the fact that 
the model is a mixture of dark matter and dark energy, we cannot exclude a priori the 
possibility that the dark energy EoS uj^ei has a phantom phase. As the observations tend 
to support, it is possible that the EoS for dark energy evolves in time and eventually 
crosses the boundary uJde = —1, at a recent epoch. In fact, as we will see below, that 
this is the case in the approach considered in this work. 



3. The Model 



Let us consider the Friedmann equation, 

H = -^P, (6) 

where H = a/ a is the Hubble parameter and p is the total density. We assume that 
this density is decomposed as p = pde + Pdm, where pde and pdm are the densities of 
dark energy and dark matter respectively. Thus for cold dark matter, the conservation 
equation becomes 

da \ f , \ /n\ 

dpde + 3 [pde + Pde) = - I "Pdm + 3 Pdm I • (7) 

The problem now is to have a relationship between dark energy and dark matter in 
order to solve ([7]). To do it, in this work we will assume a particular and simple ansatz 
which gives results in agreement with the observations, as follows, 

pie = ^Pdm, (8) 

where A is a constant. We will show that this is consistent with the usual behavior of dark 
energy, suppressed at early times, and then increasing and triggering the acceleration 
at late times. We should note that equation ([7]) has also been used in connection with 
the generalized Chaplygin gas in holographic models |16]. In contrast with our proposal 
(IH]), in these kind of models a holographic dark energy density is assumed to be able 
to relate it with the dark matter content. Under ansatz ([8]), we find the normalized 
densities 

O Pd<^ ^ JO Pd-m. Pde /„n 

^he = — = and ildm = = t— , (9) 



which clearly satisfy the fiat universe constraint Qdm + ^de = 1- If we impose the 
restriction pde > 0, then we have from 



A 

Pde - 2 



(10) 
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4. Dark energy equation of state 



Under the assumption of cold dark matter, we can identify the pressures p 
i^dePde- Therefore, 

up = UJdePde, 

and the dark energy EoS is, 

p OJ 
Ude = UJ = — . 

Pde i 'de 

After using ([2]), ([9]) and ffTOj) in the last equation, we obtain 

l+a 



^de — TT 



As 
2 



Pde = 



'111 



(12) 



(13) 



Substituting 



in this equation, we get the dark energy EoS 



2X 



l + \ 1 + 



4poXi- 
Aa3 



where 



X = l + A, [a3(i+") - 1] > 0. 



(14) 



(15) 



We can see from (1131) that Ude decreases as a increases. Moreover, for enough large 
values of a, we find in general that Ude < — 1. If we consider, the today value a = 1, we 
get 

A.. 



^deO — ~- 



1 + 



4po 
A 



(16) 



which, by considering the most probable values of As given in [15] , we show below that it 
fulfills as well the phantom divide condition Udeo < — 1- Notice that once the boundary 
LJde = — 1 is crossed, dark energy behaves as phantom-like for all the rest of the evolution 
of the Universe and it never returns to the quintessence region Ude > — 1. Moreover, in 
the limit a ~ 0, the dark energy EoS consistently approaches Ude ~ 0. 



4-1- Scaling Solution 

Scaling solutions are interesting as far as they could solve the coincidence problem 
In the scaling regime, the ratio Pdm/Pde is a non-zero constant. Thus, dark energy and 
dark matter remain of the same order. We will see that our model has an asymptotic 
scaling region for a large. In order to see it, we compute the dark energy density by 
means of (HD and 



_ A 

Pde — 2 



Aa3 



:i7) 
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Thus, taking into account ([HD we get 



Pdn 



X 
4 



1 + 



4po^ 



1 

l + ct 



Aa3 



- 1 



These are solutions of equation ([7]), taking into account f|T^ . Indeed, a simple 

calculation reproduces the GCG density (jl]). Substituting (jl]) into we get the dark 
energy normalized density in terms of the scale factor 

^ . , (19) 



de 



1 + 



Aa3 



as well as the dark matter normalized density 



1 + 



(20) 



1 + ^ 



+ 1 



Then, the ratio is obtained 



pdn 



Pde 




^PoX 



1 

l + a 



Aa3 



(21) 



which clearly differs from the non scaling models ACDM and GCG, with pdm/ Pde 
and Pdm/ Pde oc a~^''^~^"^ (see [H]) respectively. For both models cude = — 1- 
In our case, we find that this ratio is today 



PdmO 
PdeO 



1 + 

A 



(22) 



and for large values of a, we have a scaling solution 



Pde 



1 + 



A 



(23) 



We can see that (122|) and (!23l) are practically the same for the allowed values of ~ 1. 
Therefore, dark energy and dark matter are of the same order today and in the future, 
which alleviates the coincidence problem. Further, the usual "small" time interval that 
takes the transition between dark matter and dark energy (present for instance, in the 
ACDM model), in our case stretches to infinity, in such a way that the dark energy and 
the dark matter densities tend asymptotically to constant values 

-I 2 



A 

pde ^ 2 



1 + 



4po^ 



1 

l + a 



A 



- 1 



Pdv 



A 
4 



1 + 



4po^, 



l + a 
S 



A 



1 



(24) 



Let us now consider the coupled equation for the interacting model ([7]) in the scaling 



regime, when pdm and pde are practically constant. Then, dp, 



'dm 



and dpde — 0, 
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therefore pde — —{Pdm + Pde)- Thus, considering the EoS pde = ^dePdei we get in general 
for this model 



Hence, 



(25) 

Pde 



_(1 + ^)<_1. (26) 



Thus, in the scaling region dark energy EoS is phantom-like. Because the ratio pdm/ pde 

in the future is practically the same as today, we can put pdmo/Pdeo = ^dmo/^deo — 

0.3/0.7 ~ 0.43 in (1261) . Then we get the lower bound Ude — —1.43. In this limit, the 

1 

total density of the GCG model (jl]) is p ^ Ag^" po, and the Universe enters in the 

1 

de-Sitter phase a{t) ~ exp{Hft), where Hj ~ (SnGAs^" pq)/3. Therefore, the Universe 
expands accelerated forever. Now, we can analyze the coupling between dark energy 
and dark matter in more detail. We will follow the definition of the coupling given in 

m 

Pdm + SHpdm = (Hpdm, Pde + {1 + UJde) Pde = -( H Pdm, (27) 

where C is the coupling. Making the change pdm = adpdm/da, pde = ddpde/da and using 
(ED we find 

C = -3il±^, (28) 
^ l + 2r ' ^ ' 

where we introduced the ratio r = Pdm/pde- We note that both fluids always interact 

(C 7^ 0), only when Ude passes through —1/2 they uncouple. Also notice that the coupling 

evolves from a negative value through a positive value, once the condition Ude < —1/2 

is satisfied. Eventually when the approximation (125|) is valid, this coupling approaches 

to a non vanishing constant C ^ 3 for large values of a. This is consistent with the 

fact that the fluid enters in the scaling regime. It is interesting to compare our model 

with the work in [13] where Ude is considered constant. In that case the coupling (' (we 

denote with primes the quantities in that work) results in 

. -(!i±M^, (29) 

1 + r 

which for the GCG model with u'^^ = —1, one has rj = 3(1 + a) and r' = p'dml P'de — 
^dmo/^deoO'~'^ ■ Now, usiug fl27|) for our case and a similar set of equations for the model 
in [13] we find the relation between ( and (' 

C = ^^C-S — il+U^de), (30) 
Pdm dPde Pdm 

which shows that these couplings are different. Note, however, that both couplings 
coincide if Ude = — 1 and the densities with primes and those without primes coincide 
as well. 
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Figure 1. Dark energy EoS vs. loga. The solid line corresponds to As — 0.79 and 
a = 0.999. In this case dark energy EoS has a value iVdeo — —1-13 today. The dashed 
line corresponds to As = 0.936 and a — 3.75 and ujdeo — —1.34 today. In both cases 
the phantom divide is crossed recently. 

4-2. Numerical Results 

For practical purposes and in order to be specific, we shall take numerical values for the 
parameters involved in the model. Taking the present values for the fractional densities 
as fldeo = 0.7 and fldmo = 0.3, we obtain from ([9]) A = {fl'^^Q/fldmo)Po — 1.63po- The 
dark energy EoS (ITB]) is today approximately Udeo — — 1.43 A^, independent of a, and is 
phantom-like for 1 > Ag > 0.7. This falls quite well in the region of confidence given by 
the constraints on the observations for a flat Universe [I5] . We will take the best fitted 
values in the two parametrizations: a = .999, Ag = 0.79 for the range < a < 1, and 
a = 3.75, As = 0.936 for the range a > 1 [15]. In Figure 1, we show the phantom-like 
behavior of ujde- We see that it crosses the phantom divide at a ^ 0.9 for Ag = 0.79 
and has a value uJde — —1.13 today. In the case Ag = 0.936, it crosses at a ~ 0.95 and 
has a value uJde — —1.34 today. In Figure 2, we plot the quotient of the dark matter 
over dark energy densities for the same values as in Fig. 1. We see that it approaches 
to a non vanishing constant limit as the Universe expands, and it remains constant in 
this region. Further, from (l23ll we have in this region, pde — 2.55 pdm for Ag = 0.79 and 
Pde — 2.35 Pdm for Ag = 0.936. In Figure 3, we plot the coupling for the same choice of 
the parameters. As we can see, the coupling is constrained to C < 3. Finally, in Figure 
4 we plot Qdm and Qde- In the past they behave as in the ACDM model, but in the 
future for Ag = 0.79 they approach to the limits Qdm — 0.282 and fide — 0.718 whereas 
for Ag = 0.936 they approach to Qdm — 0.298 and Qde — 0.702. These are almost the 
same values as the ones at present. 
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Figure 2. Dark matter density over darlc energy density r vs. logo . The solid line 
corresponds to — 0.79 and a — 0.999. The dashed line corresponds to Ag = 0.936 
and a = 3.75. 
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Figure 3. Behavior of the coupling C vs. loga. The solid line corresponds to a = 0.999 
and As = 0.79 and the dashed line corresponds to A^ = 0.936 and a — 3.75. 

5. Conclusions 

In this work we have studied the generalized Chaplygin gas as an unified self-interacting 
fiuid. By assuming a simple relation between dark matter and dark energy densities ([H]) 
we obtain a dark energy EoS which exhibits a phantom-like behavior. In the regime of 
the most probable values of and a [15] , this identification leads to a dark energy EoS 
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log(fl) 

Figure 4. Normalized dark matter ^Idrn (decaying curves) and dark energy il.de 
(growing curves) vs. loga. The solid lines correspond to As = 0.79 and a = 0.999. 
The dashed lines correspond to Ag = 0.936 and a = 3.75. We see that they behave 
like the ACDM model in the past. 

that crosses the phantom regime = — 1 in the recent time. In fact, the numerical 
values Udeo — —1.13 (for Ag = 0.79) and Udeo — —1.34 (for Ag = 0.936) are consistent 
with the observations. Furthermore, the general behavior of the dark energy EoS in 
our model, is in complete agreement with the best fit in [iTj for Qdmo = 0.3. Once the 
dark energy EoS crosses the phantom barrier, it never returns to quintessence, giving a 
phantom-like sector for all the future evolution of the Universe. The densities smoothly 
decay for all values of the scale factor, thus avoiding a Big Rip. We find also, that in 
contrast with the ACDM and GCG with Ude = —1 models |13j, our model exhibits a 
scaling solution when a is large, as can be seen from the numerical results. This regime 
is achieved once the dark energy EoS has crossed the phantom divide, and remains in 
this way for all the rest of the evolution. We show also that the coupling between dark 
energy and dark matter tends to a positive constant in this region and the Universe 
tends to be de-Sitter. This type of solution could provide a clue to solve the coincidence 
problem [Hj. We also find that the present model with time varying EoS and with 
the assumptions made, could be considered at least as a phenomenological model which 
could encode the interaction conjectured in [13], inferred from the Abell cluster A586. 
The relation between both interacting models is given through (l30l) . It will be of interest 
to find scalar fields and their corresponding potentials that could reproduce our results. 
This will be studied in a future work. 
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